Distribution of extremes in the fluctuations of two-dimensional equilibrium interfaces 



Deok-Sun Lee 

Theoretische Physik, Universitdt des Saarlandes, 66041 Saarbrucken, Germany 
(Dated: February 2, 2008) 

We investigate the statistics of the maximal fluctuation of two-dimensional Gaussian interfaces. Its re- 
lation to the entropic repulsion between rigid walls and a confined interface is used to derive the average 
maximal fluctuation (m) ~ \j2/{nK) h\N and the asymptotic behavior of the whole distribution P(m) ~ 
pjl e - (const) N-e m -\f2nKm f Qr m fj n i te w j m jy2 an£ j g ^ m t er f ace s j ze tension, respectively. The stan- 
dardized form of P(m) does not depend on N or K, but shows a good agreement with Gumbefs first asymptote 
distribution with a particular non-integer parameter. The effects of the correlations among individual fluctua- 
tions on the extreme value statistics are discussed in our findings. 
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A common nature of floods, stock market crashes, and the 
Internet failures is that they hardly occur but, once they do, 
have significant consequences for the corresponding systems. 
These rare but fatal events correspond to the appearances of 
extreme values in the fluctuating variables such as the daily 
discharge of a river, the stock index, and the router load, and 
therefore it has been in great demand to be able to estimate 
and predict the typical magnitude of the largest values and the 
probability of a given extreme value 01 • The extreme value 
theory has been of great importance also in various physics 
contexts since the physics of complex systems is often gov- 
erned by the extreme value statistics |2]. Recently, Bramwell 
and co-workers discovered a universal distribution for the or- 
der parameters of a class of model systems for ferromagnet, 
confined turbulent flow, avalanche, and granular media |3|,|4]]. 
The universal distribution turns out to be consistent with Gum- 
bel's first asymptote distribution that is one of the well-known 
probability distributions, to be precise, of the nth largest value 
among 0\[ independent (uncorrelated) random variables |1]. 
Interestingly, the parameter n takes a non-integer value and 
this has been considered as related to strong correlations in 
those systems |4j]. While the reason for this universality is 
still debated 1 5, 6], it renewed the interest in the extreme value 
statistics of correlated variables 0001. 

In this Letter we study the statistics of the maximal fluc- 
tuation of two-dimensional (2D) equilibrium interfaces. The 
Gaussian ensemble is considered, where an interface con- 
figuration {(j)} = {(j)(r)|r = (x,y),—L/2 < x,y < L/2] un- 
der the periodic boundary condition has the probability 
P[{ty}] = e^oIWl with the Gaussian Hamiltonian #b[{0}] = 
/ d 2 r(K/2)(V§) 2 and K the interface tension. The spatial 
average (]) = L~ 2 J d 2 r§(r) is set to be zero. This Gaus- 
sian model has been used to study rough (self-affine) in- 
terfaces ubiquitous in nature [ 10], the interacting spin sys- 
tem at low temperatures through the 2D XY model with 
(j) representing spin orientations |4], and so on. We de- 
fine the maximal fluctuation m for a given interface {(])} as 
m = maxj! G /f{|(|)(r)|}, where/? = {(an x ,an y )\n x ,n y = —N/2 + 
l,—N/2 + 2,...,N/2} with a the lattice constant and N = 
L/a. It is an extreme value among the individual fluctua- 
tions, (j)'s, that are correlated via as ((j)(r)(|)(r')) ~ln(L/|r— 



V\) fioll . and we study the statistics of m. The effects of such 
correlations have been studied by a renormalization group ap- 
proach in the context of glass transition 1 11], but the nature 
of the statistics of m remains to be understood. The maximal 
fluctuation has its own technological significance as well, for 
instance, in the onset of the breakdown of corroded surfaces 
and the occurrence of a short circuit by the metal surface of 
one electrode reaching the opposite one fl2lfl3ll . 

The exact distributions of the maximal height fl2l [l3ll and 
the width square 1 14] have been obtained for the ID Gaussian 
interface. For two dimensions, however, the distribution of 
the maximal fluctuation is not known while the width-square 
distribution is well understood 1 15]. Of our particular interest 
are the functional form of the maximal fluctuation distribu- 
tion, P(m), and its first few moments. We use the relation of 
P(m) to the entropic repulsion between an interface and rigid 
walls to find (m) = J dm mP{m) and the asymptotic behavior 
of P(m), which is shown to resemble Gumbel's first asymp- 
tote distribution. The standardized distribution of the maximal 
fluctuation does not depend on any system parameter, and fur- 
thermore, Gumbel's first asymptote distribution with a spe- 
cific parameter fits excellently the standardized distribution of 
m. These results are compared with the statistics that would 
appear without the correlations among (|)(r)'s to illuminate the 
correlation effect on the extreme value statistics. 

When an interface is confined between two rigid walls, the 
free energy is increased due to the entropic repulsion against 
the walls. The maximal fluctuation distribution P(m) of a free 
interface is related to this free energy increase with the walls 
at (]) = m and (j) = — m, Af (m), as 
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where % m [{0}] is 1 f° r max ?e/f{| ( K' i )|} — m an d otherwise, 
and thus XmlW] = I\reR®( m — with 9(jc) = 1 for x > 

and otherwise. We first derive Af (m) and obtain P(m) 
using Eq. Q. 

For the soft walls restricting the width square W2 = §(r) 2 to 
be less than x 2 , the free energy increase is given as Af (x) ~ 
The entropy reduction is inversely 
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proportional to the projected area of the interface segment 
I 2 that is so large that the average fluctuation over the seg- 
ment, ~ Vln I, is comparable to x. For the rigid-wall prob- 
lem, one may expect the same result if the average fluctu- 
ation and the maximal fluctuation scale as functions of the 
segment area in the same way, which is true for ID inter- 
face but not for the 2D one. Bricmont et al. have 
shown, by means of a rigorous proof as well as an heuristic 
argument, that Aj (x) ~ e - (const)* f or confined Gaussian 
interfaces in the thermodynamic limit N — *■ °° 1 18]. Reviewing 
briefly the heuristic argument introduced therein, we compute 
Af (x) keeping the lateral size N large but finite, which leads 
us to find out (m) and two distinctive behaviors of P(m) for 
m -C (m) and m ^S> (m) , respectively. 

One of the most prominent effects of confinement on an in- 
terface is the reduction of the correlation length; the height 
fluctuations are hardly correlated across the interface-wall 
collisions. The correlation length can be obtained in a self 
consistent way by assuming that the effective Hamiltonian 
#"eff(je) satisfying Tre -5 ^*) /Tre - * = e~ A!r W with Aj (x) 
in Eq. Q takes the form of a massive Gaussian Hamiltonian, 
X eS {x) = Jd 2 r[(K/2)(y§) 2 + (ju(x) 2 /2)$(7) z }. The effective 
mass fi(x) is the inverse of the correlation length and obtained 
self-consistently. Differentiating Eq. Q with respect to m = x 
and using !H e fs{x), one obtains an approximate expression for 
the derivative of the free energy increase as 
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where e ^°%x is approximated by e H &\ x ) and (A) f( s t x \ = 
TrAe~*" ef M/ Tre ~*° eE ^- Also we used the relation that 
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For Eq. 0, we considered the Fourier components fa = 

L- 2 Jd 2 re- ( ^^(r) for q = (n/L)(2n x - l,2n y - 1), which 
are decoupled in ^{ B g(x), allowing another expression for 

— (d/dx)Aj (x) as 

3 , , d f Tre-^ffW\ 
- —A 7 (*) ~ — In — 
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Equations (|2j-@ give a self-consistent equation for n(x), 
which can be solved numerically with the condition n{°°) = 0. 
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FIG. 1: Values of -3Ay (jc) /dx with K = 0.92 calculated by substi- 
tuting P(m) from the Monte Carlo simulation to Eq. [N = 32(D), 
64(o), 128(A)] and by the self-consistent solution to Eqs. 0-0 
[N = 32 (dotted line), 64 (dotted-dashed), 128 (dashed)]. 
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FIG. 2: (a) Plots of ln(ra) and In W versus ln(lnA') for 2D free Gaus- 
sian interfaces with K = 0.50, 0.92, and 1.50. The upper solid line 
has the slope 1 while the lower dashed line has 1 /2. The dotted lines 
connecting data points are guide for the eye. (b) Semilog plot of (m) 
versus N for the same values of K as in (a). The lines have the slopes 
y/2/{nK), respectively. 



A comparison of the values of — dAf (x)/dx from this equa- 
tion and from a Monte Carlo simulation for free interfaces 
giving P(m) and in turn, — dAf (x) /dx through Eq. Q, can 
be a test for the validity of the adopted approximation. As 
shown in Fig. both are in excellent agreement except for 
a slight constant deviation that is presumably due to higher- 
order (in (]) 2 ) contributions to the singular factor Q(x — \§\) = 



limr, 



-(♦A)' 



1 18]. This justifies the use of ?4ff(x) for the 
study of the maximal fluctuation distribution of a free inter- 
face as well as the free energy increase of a confined interface. 

The large-* behavior of the effective mass fi(x) can be 
derived analytically. From Eq. 0, W(x) 2 is represented 
as (2nK)- 1 ln\p Q /ij(x)} for no/fi(x) <C N and (2nK)- 1 InN 
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for jUo//i(jc) ^S> N, where /jo = \[K%ja. Solving Eqs. (|2j 
and @ with these values of W(x) 2 , one finds that, for 
large x, /j(x)//jq ~ e -^J nK l 2x forl«i«i f and fi(x)/fio ~ 
g-nKf/pbxfr) foix > ^ with ^ _ ^/2/{%K) In AT. It is natural 
to think that the interaction between the walls and the interface 
is negligible and that the correlation length diverges when the 
distance to the wall x is much larger than (m). Therefore, we 
expect the characteristic distance x c is equal to (m), 



— In N. 
%K 
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This relation is in agreement with the result in Ref. full and 
is also verified by our Monte Carlo simulation data for (m) 
as a function of N plotted in Fig. [2] The scaling of (m) is 
distinguished from that of the roug hnessW- VlrWV[Fig.|2l. 
This is contrasted to the ID Gaussian interface, where both the 
maximal fluctuation and the roughness scale as y/N Jl2T.[l3ll . 

Using the obtained functional form of /j(x), one can also ob- 
tain the free energy increase Af (x) and in turn, the maximal 
fluctuation distribution P(m) through Eq. Q: 




FIG. 3: Skewness of the Gumbel distribution Yg (' J ) = / dzz 3 g(z;n). 
It is shown that Yg(") with n = ;;,„ ~ 2.6(2) is equal to the skewness 
of P(m), y m , which is 0.68(2) independent of the system size N or 
tension K as shown in the inset. The skewness of the width-square 
distribution, y Kl , is about 0.89, which is reproduced by g(z',n) with 
n = n w ~ 1.58H. 
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The asymptotic behavior of P(m) for m <C (m) in Eq. © is 
consistent with Gumbel's first asymptote distribution that rep- 
resents the limiting distribution of the n th largest value, m n , 

among 9\[ — > °° independent random variables following an 

l"T, 

exponential-type distribution 1 1]. In terms of the standardized 
variable z — (m„ — (m„))/a,„ n with (m„) and G m „ the aver- 
age and the standard deviation of m„, Gumbel's first asymp- 
tote distribution takes the standardized form g(z;n) given 
by OH g(z\n) = CO e -nle- b ^+b{z+ S )] f where b = ^y^), 
s = (Inn — \y{n))/b, and CO = n n b/T{n) with T(x) the Gamma 
function and \|/(jc) = 8lnr(x) /dx 01 . 

It is g(z; n) with a non-integer value of n that fits the univer- 
sal distribution in a class of correlated systems The 
2D Gaussian interface belongs to this class since the stan- 
dardized distribution of the width square Q W2 (z) = W2 P(\V2 = 
(W2) +zG W2 ) with (w2) and a W2 the average and the standard 
deviation of W2, respectively, is fitted by g(z;n — ~ 1.58) 
without regard to the system size N or tension K |3l|4|]. To 
check whether such a unique distribution also exists for the 
maximal fluctuation, we computed the skewness of P(m), 
y m = ((in — (m)) 3 ) /af n with a,„ the standard deviation of m. 
According to our data shown in Fig. [3] y m hardly varies with 
N or K, but stays at 0.68(2), which leads us to expect a sim- 
ilar universality in the maximal fluctuation statistics to in the 
width-square one. That particular value of the skewness y m 
is reproduced in Gumbel's first asymptote distribution with 
n = n m ~ 2.6 as shown in Fig. [3] and we compare the stan- 
dardized distribution Q m (z) = O m P(m = (m) +zO m ) for dif- 
ferent values of N and K with Gumbel's first asymptote dis- 
tribution g(z;n m ) in Fig.|4] One can identify the agreement of 
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FIG. 4: Standardized distribution Q m {z) for N = 32, 64, 128, and 
K = 0.5, 0.92, 1 .5. The data points in the range z < 5 fall on the Gum- 
bel distribution g(z;n m ) with n m ~ 2.6 represented by the solid line. 
Inset: The standard deviation o" m of the maximal fluctuation is fitted 
by 0.29(1)Jl -1 / 2 (solid line) without a significant A'-dependence. 



Qm(z)'s and g(z;n m ) as well as the data collapse of Q m (z)'s 
for different APs and K's. In addition, this agreement en- 
ables one to predict that o,„ is given by o„, = b/y/2%K with 
b = yW^fhn) — 0.69(3), which is confirmed numerically in 
Fig.E 

Now let us look into the effects of the correlations among 
<j)(r)'s on (m) in Eq. I0 and on P(m) in Eq. (|6j. For compari- 
son, we consider the statistics of m that would emerge with no 
correlation for a Gaussian distribution of 5\£ individual fluctu- 
ations p(§) =e~§ <( 2W I /V2nW 2 . In this case we could apply 
the following relation to obtain (m): JZ^Ldtypfy) ~ — 
Vj/0t (H. This relation yields (m) ~ W^jMpT/W) O. The 
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yln3\£ -scaling of (m) is thus generic when the roughness W 
is finite. A different scaling of (m) shows up with a diverging 
roughness. For instance, (m) ~ ^/N\nN in ID Gaussian in- 
terfaces where 9\[ = N and W ~ \/N and (m) ^h\N in 2D 
ones where 9£ = N 2 and W ~ VlniV. While the y/Nh.N 
behavior deviates from the true behavior of (m) in the ID 
Gaussian interface lll2l[l3ll . the ln^V behavior in two dimen- 
sions is consistent with Eq. £J. Without the correlations, 
the distribution of the maximal fluctuation could also be ob- 
tainedby jl| P{m) ~25\£ (J m m d§p(§))^ p(m). This gives 

P(m) ~ 91 e -(const)^-<" 2 /( 2,v2 }-m 2 /(2W 2 ) ^ p ^ for m ^ ^ 

but is not capable of reproducing the behavior of P(m) for 
m finite in Eq. 0, which demonstrates the essential role of 
the correlations. The extreme statistics without correlations 
has been found in small-world networks where the correla- 
tion length is finite due to random links fl9il . The correla- 
tion length becomes finite also in the presence of external 
field 12011 . P(m) is more skewed, that is, the skewness is larger, 
with larger external field |21]. On the contrary, the width- 
square distribution recovers its Gaussian form as it loses the 
correlations 1 20] . 

Our results demonstrate the availability of Gumbel's first 
asymptote distribution with a non-integer parameter n for the 
maximal fluctuation distribution of the 2D Gaussian inter- 
face. The value of n can be different according to the type 
of the correlations: It takes a value close to k/2 in the case 
of the maximal avalanche distribution in the Sneppen depin- 
ning model [8], different from 2.6(2) found for the maximal 
fluctuation distribution in this work. We remark that the value 
of n m is larger than n W2 , that is, y W2 > J m , implying that the 
width-square distribution is more non-Gaussian than the max- 
imal fluctuation distribution. 

In summary, we investigated the statistics of the maximal 
fluctuation in the 2D Gaussian interface. Because of the cor- 
relations among individual fluctuations, the average maximal 
fluctuation as well as the whole distribution exhibit distinct 
behaviors from those of uncorrected fluctuations. We also 
identified the appropriateness of Gumbel's first asymptote dis- 
tribution with a non-integer parameter for the maximal fluctu- 
ation statistics, which suggests its applicability to the extreme 
value statistics in a class of correlated systems with the pa- 



rameter depending on the correlations. 
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